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Abstract. We present a covariant approach to describe the low-lying scalar, pseudoscalar, vector 
and axialvector mesons as quark-antiquark bound states. This approach is based on an effective 
interaction modeling of the non-perturbative structure of the gluon propagator that enters the quark 
Schwinger-Dyson and meson Bethe-Salpeter equations. We extract the meson masses and compute 
the pion and kaon decay constants. We obtain a quantitatively correct description for pions, kaons 
and vector mesons while the calculated spectra of scalar and axialvector mesons suggest that their 
structure is more complex than being quark-antiquark bound states. 

INTRODUCTION 

Recently, the scalar mesons have attracted a lot of interest as the reanalysis of the 
pseudoscalar meson scattering data indicated the existence of a flavor SU(3) nonet in 
this channel 1 1]. It is therefore desirable to gain deeper understanding of the constituent 
structure of the scalar mesons together with a comprehensive description of the meson 
states in the other spin-parity channels. The ultimate goal would be to understand all 
low-lying meson states and resonances as non-perturbative bound states in Quantum 
Chromo Dynamics (QCD). 

A relativistic framework for analyzing mesons as composite objects is provided by the 
Bethe-Salpeter equations that extract poles in the quark-antiquark scattering kerneljsl]. 
The attraction needed to bind quarks and antiquarks emerges from dressed multiple 
gluon exchange. Thus the essential ingredients to these equations are the quark and 
gluon propagators as well as the quark-gluon vertex. In addition, these n-point Green's 
functions are related by their Schwinger-Dyson equations which are part of an infi- 
nite tower of non-linear integral equations. There has been some progress in the under- 
standing of the infrared behavior of the gluon propagator from recent Yang-Mills lattice 
measurements |2] as well as from studies of the coupled system of gluon and ghost 
Schwinger-Dyson equations 10, 0]- More recently, the coupled system of gluon, ghost 
and quark propagators functions has been studied within certain truncation schemes of 
the Schwinger-Dyson equations and solutions have been obtained for various ansatze 
for the ghost-gluon vertices [5]. Nevertheless, for phenomenological applications the 
frequently adopted strategy is to model the gluon propagator as well as the quark-gluon 
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vertex and consistently derive the quark propagator from its Schwinger-Dyson equation. 

These types of calculations have a long history, for reviews see refs. [3, 6]. Early 
versions adopted pointlike gluon propagators in coordinate space that eventually lead 
to Nambu-Jona-Lasinio (NJL) type models iQSI], pointlike propagators in momentum 
space were also considered [ja]- These models are particularly simple because either 
solving the Schwinger-Dyson equation yields a free quark propagator or the Bethe- 
Salpeter integral equations reduce to algebraic equations. The main target particularly of 
the NJL-model studies have been the pseudoscalar mesons. It turned out that they can be 
adequately described once the important feature of dynamical chiral symmetry breaking 
is incorporated, i.e. the interaction is strong enough so that the resulting quark propagator 
develops a non-zero constituent quark mass. Then the pseudoscalar mesons can be 
understood as the would-be Goldstone bosons of chiral symmetry breaking. However, 
these NJL-type models do not reflect the confinement property of QCD and thus binding 
can only be achieved kinematically, i.e. meson states with masses larger than twice the 
constituent quark mass cannot be described consistently. For that reason, model gluon 
propagators have been developed that yield quark propagators without poles for real 
momenta as an attempt to include the confinement phenomena fl^ [Tl|]. Again these 
studies focused on pseudoscalar mesons [12, 13] while a comprehensive investigation 
for the scalar, pseudoscalar, vector and axialvector mesons has not been carried out so 
far. Other studies ifToLITll] made contact with perturbative QCD by considering a model 
gluon propagator that matches the pertinent anomalous dimension. This contribution 
has negligible effect on the meson properties, but its inclusion makes cumbersome the 
extraction of the solutions to the Schwinger-Dyson equations for the large time-like 
momenta that enter the Bethe-Salpeter equations. Such large time-like momenta need 
to be considered for mesons other than the pseudoscalars. Although this is interesting 
we regard it an unnecessary technical complication because we do not want to compute 
properties of mesons revealed only at high momentum transfer. Rather we want to 
establish a model as simple as possible that we consider a pertinent starting point to 
study the structure and properties of low-lying mesons in a fully relativistic framework. 
Our model interaction is parameterized in form of a non-trivial gluon propagator that 
contains sufficient strength to cause dynamical chiral symmetry breaking. For technical 
reasons it turns out that a GauBian shape function for the propagator in momentum 
space is most suitable. Essentially we consider this model propagator as an effective 
interaction that relativistically describes the binding of quarks and antiquarks to mesons. 
Furthermore, we take the quark-gluon vertex function to be the tree level one since this 
procedure provides a framework that is consistent with chiral svmnietry when the ladder 
approximation for the Bethe-Salpeter equation is employed ylla]- For approaches going 
beyond ladder approximation see e.g. ref. 

This talk, which is mainly based on ref. [J^, is organized as follows: First we will 
introduce the effective interaction and solve the Schwinger-Dyson equation for the 
quarks. We will put particular emphasis on the analytic continuation of the resulting 
quark propagator to time-like momenta that enter the Bethe-Salpeter equations. We will 
discuss the structure of the Bethe-Salpeter equations and then present solutions. Finally 
we will conclude and suggest a possible extension of the current approach in particular 
with regard to the possibility that the scalar meson might have to be considered as two- 
quark - two-antiquark bound states lIlLfl^. 



THE QUARK SCHWINGER-DYSON EQUATION 



We take a GauBian form for dressing the model gluon propagator and write, 

g^Gf,{q) = An^D5"\,{q)^ exp [-^ (1) 

where /i, v are Lorentz indices, tp_v{q) is the transverse momentum projector and a,b 
label color. While the coefficients in eq. ^ are chosen to make subsequent equations 
more concise, D and O) are dimensionful parameters that we will determine from fitting 
empirical data. The coefficient D sets the strength of the interaction and O) is the value at 
which the scalar function in the parameterization is maximal. Hence CO sets the interac- 
tion scale. The dressed gluon propagator dU) is supposed to represent a sensible hadron 
model and hence one can envisage that CO will have a value of several hundred MeV. 

We interpret the effective interaction ([T]) as the propagator (in Landau gauge) of a 
gluon that gets absorbed and emitted by the quarks that eventually get bound to form 
mesons. To completely define the interaction, we need to parameterize the quark-gluon 
coupling. To establish chiral symmetry we apply the rainbow-ladder approximation to 
the system of Schwinger-Dyson and Bethe-Salpeter equations. Then the quark-gluon 
coupling is given by the tree level interaction vertex, igf^^, where A" is a Gell-Mann 
matrix acting in color space. Note, that we have already included the coupling constant g 
in the definition of the effective interaction dl]). 

The Schwinger-Dyson equation for the (inverse) quark propagator becomes 

S-\p)=ip + mo + J ^,ynS{k)y^g'^^Gf,{k-p) (2) 

where niQ is the current mass of the considered quark. This contribution represents the 
only explicit distinction between quarks of different flavors. Of course, its effects will 
implicitly propagate through the whole calculation. However, for notational simplicity 
we will continue to suppress flavor labels. A suitable parameterization of the quark 
propagator is inspired by the form of a free fermion propagator 
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In solving the Schwinger-Dyson equation Q we have to find the scalar functions A (p^) 
and B{p^). It is also very instructive to define a mass function via 

M(p2) =B{p^)/A{p^). 

In particular M{p^ = 0) plays the role of a constituent quark mass and a large value 
(^ mo) thereof signals dynamical chiral symmetry breaking. 

We work in Euclidean space with Hermitian Dirac matrices that obey {y^x, 7v} = 25^v 
and 75 = —71727374. Inserting the effective interaction ([T]) and performing the standard 
trace algebra, we then deduce the following coupled equations for the propagator func- 
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FIGURE 1. Quark propagator functions (as a function of the momentum squared). The parameters are 
CO = O.SGeV, D = 16.0GeV"^. Left panel: otq = 0, right panel: mo = 0.115 GeV. All units are GeV. 
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where the four dimensional integral measure has been expanded such that x = and 
y = k^. Furthermore /„ are modified Bessel functions. Note, that as a particular feature of 
the GauBian dressing function ([T]) it has been feasible to compute the angular integrals 
analytically. 

In a first step we solve eqs. Q for spacelike momenta, i.e. for real positive x. Then 
we observe that the integrals on the RHS of these equations only involve the propagator 
functions at real arguments y and we can use them to numerically compute the propaga- 
tor functions for arbitrary complex x. At first sight, it appears that A(x) and B{x) could 
not be consistently continued because the cut along the negative x-axis (associated with 
^/x) would yield different results when continuing in the upper or the lower half-plane 
and it would be impossible to resolve the ambiguity in ^ ^ when continuing 

X — i> — Fortunately, this is not an obstacle because the modified Bessel functions I\ 
and I2 are respectively odd and even functions of their arguments. Thus we are free to 
choose either of the two signs above. For definiteness we work in the upper half-plane 
with y/x —>■ i\/\ along the negative half-line. 

In Fig. [U we show the quark propagator functions A(p^), and M{p^) — 

B{p^)/A{p^) along the positive, real spacelike axis, p^ > 0. This is the numerical so- 
lution to the coupled equations which we emphasize is the basis for the quark so- 



lutions for complex momenta. The solution clearly shows that dynamical chiral sym- 
metry breaking is occurring: the mass function M{p^) attains a sizable non-zero value, 
even in the case that the bare quark mass mo is zero. This phenomenon is an example 
of genuinely non-perturbative behavior as dynamical mass generation cannot occur at 
any order in perturbation theory. Recall that the effective interaction ([T]) that enters the 
Schwinger-Dyson equations does not contain the perturbative UV behavior, rather it 
has an exponential damping at high momenta. This is manifested in the quark propaga- 
tor functions as a sharp transition from the low momentum behavior to the bare values 
in the high momentum region. This transition occurs at about IGeV. 



THE BETHE-SALPETER EQUATION 

Having obtained the quark propagators in the complex plane from the Schwinger-Dyson 
equations we have collected all ingredients for the Bethe-Salpeter integral equations. 
They will ultimately yield the quark meson vertex functions, F, that describe mesons 
as bound quark-antiquark pairs. These vertex functions for bound states are determined 
from the condition that the four quark scattering amplitude develops a pole at = —M^, 
where M is the mass of the bound state meson. 

The vertex functions resulting from the Bethe-Salpeter equation are characterized 
by three momenta out of which only two are linearly independent due to momentum 
conservation at the vertex. If we denote the meson momentum P and the momentum of 
the incoming quark p + ^P then the momentum of the outgoing quark (= incoming 
antiquark) is p + — l)P. This suggests to label the vertex functions by p and P: 
r{p, P) . We have introduced the arbitrary momentum partition parameter £, G [0, 1] . Due 
to strict relativistic covariance the results for physical observables do not depend on_£ . 
We have studied the (^-dependence of our numerical results exhaustively in ref. [ 15ll . 
The confirmed (^-independence represents an a posteriori validation for the relativistic 
covariance of our computations. 

We now turn to the main target of our studies, the Bethe-Salpeter integral equations 
for the vertex function Y{p,P) in ladder approximation I3..6. .17.1 : 

r(p;/') = -^/^ [yvS{k+^p)nk-p)s{k+{^-i)p)y^]g'G^^{k-p). (5) 

Here we have factorized the color factors in the effective interaction, G'^y{q) = 

8"^Gpiv{q) and performed the corresponding trace. The flavor content of the meson 
is not made explicit in eq. Q as we have suppressed the flavor labels in the quark 
propagators. It is understood that the two propagators in eq. ^ are taken such as to 
account for the flavor quantum numbers of the considered meson. In the model that we 
will consider, the up and down quarks will be assumed to have equal current masses 
(mo in eq. Q) and thus also identical propagator functions A{x) and B{x). For the 
light quarks, which should give rise to the familiar SU (3) / nonet, we are thus left with 
three representatives of each of the multiplets that are distinguished by their isospin 
number, / = 0, 5, 1 . We must also specify the meson angular momentum and parity. 
This is reflected by the Dirac and Lorentz decomposition of the meson vertex functions. 



This decomposition is known in the literature and here we follow ref. ifisll . For the 
pseudoscalar channel (7^ = 0^) we take 



r(^)(p;P) = [r^^^\p;p)-ipr[''\p;P)-ipri''\p-p) - [p,/>]rP{p;P) 

The decomposition for a scalar (7^ = 0+) meson reads 
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The vector (7^ = 1 ) channel involves eight scalar functions 
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In the axialvector channel we have two modes that are distinguished by their charge 
conjugation properties [18]. For 7^*^ = 1++ the decomposition is 
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mode is decomposed as 
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In what follows we will omit the superscripts that label the spin and parity channels 
because these channels do not mix and there should hence be no confusion. 

The solution to the Bethe-Salpeter equation not only yields the meson masses but also 
the meson quark vertex functions that can be used to compute meson properties. Here 
we will focus on the pseudoscalar decay constants f^i and fx- In order to calculate these, 
we first have to normalize the vertex functions Y{p;P). The Bethe-Salpeter equation 
is a homogeneous equation, and thus needs an additional normalization condition. As 
mentioned previously, that condition is obtained from demanding the pole in the four- 
quark Green's function to be unity. For equal momentum partitioning, (i.e. for = 1/2 
only) it reads iUtI] 
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where the trace is over Dirac matrices. The conjugate vertex function T is defined as 

T{p,-p)=cr^{~p-p)c-\ (12) 

where C = —7274 is the charge conjugation matrix. The quark propagator derivatives are 
calculated by differentiating the quark Schwinger-Dyson equations Q analytically and 
then numerically integrating the corresponding expressions. 

The decay constants are finally obtained from the coupling of the axial current to the 
quark loop [.17i] 

f=^J ^Tr{r{k-P)S{k + P/2)r5PSik-P/2)}. (13) 

The primary subject of this talk is to extract the bound state masses for the various 
flavor combinations and angular momentum channels. The corresponding projection 
results in sets of coupled equations for the F, . After carrying out two of the three angular 
integrals analytically we are left with functions of the squared momenta and P^ as 
well as the angle between p and P: z = p P/ \/ p^P^. The z-dependence is analyzed by 
an expansion in Chebyshev polynomials Tj^ 

T,{p-P)=Y,iifT\{p\P^mz). (14) 

k 

Since the T]f_ form an orthonormal set, we can project the equations for the Dirac com- 
ponents onto r^{p^;P^). Finally, the ^^-integral in the Bethe-Salpeter equation ^ is 
implemented numerically as a matrix equation for the unknown F^(p^;/'^), pj being the 
discrete values of the momentum squared. The kernel, K of that matrix parametrically 
depends on the meson momentum P^. We solve that matrix equation as an eigenvalue 
problem by tuning the meson momentum to P^ = —M^, such that Det(l — ^) = 0. This 
yields the desired meson mass M. 

In the discussion of numerical results we note that there are four model parameters, 
(O.D.niu and that we first have to fit to empirical data. To this end, we initially 
choose the pseudoscalar meson observables Mji^Mk and /^r- Then one parameter remains 
unconstrained by the pseudoscalar sector alone. However, the condition that the quark 
propagator function reflects dynamical chiral symmetry breaking, i.e. that M{p^ = 0) ~ 
0.5GeV leaves only a small window for the remaining choice. All other masses and 
decay constants are subsequently model predictions. The resulting model parameter and 
the predicted kaon decay constant fK, that is unexpectedly well reproduced, are shown 
in Table [U The subsequently predicted meson masses are shown in Tables |2l|5l In all 
cases we have the inequalities Muu < Mus < Mss, where the subscript labels the flavor 
content. These relations just reflect the quark-antiquark picture that is implicit in the 
present Bethe-Salpeter approach. 

Obviously both the pseudoscalar (table [l]) and vector mesons (table ^ can be very 
well described within our model with the choice CO ~ O.SGeV. Our results agree with a 
previous analysis of the vector mesons based on an effective interaction which included 
the perturbative type term [11]. This shows that such terms do not have a large effect 
on the meson masses, at least for the pseudoscalar and vector cases. Indeed, in the 



TABLE 1. Parameter sets used and fit results for the pseudoscalar 
mesons. M-ji, f-ji and Mk are used as input, fK is predicted. All units are 



GeV. 



CO 


D 


'«« 


m.s II Mji 


h 


Mk 


h 


0.40 


45.0 


5 X 10-3 


0.120 


0.135 


0.131 


0.496 


0.164 


0.45 


25.0 


5 X 10-3 


0.120 


0.135 


0.131 


0.496 


0.163 


0.50 


16.0 5x10-3 


0.115 


0.137 


0.133 


0.492 


0.164 


experiment 1191 || 0.135 


0.131 


0.498 


0.160 



TABLE 2. Results for the vector mesons. All units are GeV. 



(0 


D 


m„ 


rris 


Mp 


Mk* 




0.40 


45.0 


5 X 10-3 


0.120 


0.748 


0.939 


1.072 


0.45 


25.0 


5 X 10-3 


0.120 


0.746 


0.936 


1.070 


0.50 


16.0 


5 X 10-3 


0.115 


0.758 


0.946 


1.078 


experiment 1191 




0.770 


0.892 


1.020 



context of low-energy meson phenomenology we conclude that the logarithmic tail, 
and its associated renormalization represent an unnecessary obfuscation. 

The situation for the scalar mesons (table |3l) is not quite that clear. To begin with the 
particle data group [19] does not provide a clear picture in this channel but only quotes 
a wide range for the mass of the lowest scalar (0.4 - 1.2 GeV). More detailed studies of 
the pseudoscalar scattering amplitudes revealed that the assignment of the scalar meson 
nonet is not at all established [1]. In particular, these mesons may not be simple quark- 
antiquark bound states but e.g. might contain sizable admixture of 2quark-2antiquark 
pairs lH^ . In that respect we might interpret our results as a quark-antiquark model 
prediction for scalar mesons. Our results suggest that such a picture is too simple for 
these mesons. One might also speculate that the adopted ladder approximation could be 
insufficient. 

For the axialvector mesons we have two channels that are distinguished by their 
charge conjugation properties, cf. tables |4l and |5l The quark-antiquark pairs that are 
bound to axialvector modes with negative charge conjugation eigenvalue tend to be 
lighter than those with the positive eigenvalue but otherwise equal quantum numbers. 
Generally we find that our predictions are lower than the assignments made by the 
particle data group lfl9ll . 

We recognize from our results that the model predictions change only slightly within 
the large range of considered model parameters. This confirms that meson static proper- 



TABLE 3. Results for the scalar mesons. The subscripts of 
M denote the flavor content. All units are GeV. 



(0 
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niu 


nis II Miiu 


Mui 




0.40 


45.0 


5 X 10-3 


0.120 


0.700 


0.917 


1.096 


0.45 


25.0 


5 X 10-3 


0.120 


0.675 


0.908 


1.099 


0.50 


16.0 


5 X 10-3 


0.115 II 0.645 


0.903 


1.113 



TABLE 4. Results for the axial-vector (7 *' = 1+^) mesons. 
The question mark indicates that the PDG did not assign the 
charge conjugation property of the respective resonance. All 
units are GeV. 



CO 


D 


mu 


m.s II M„T, 






0.40 


45.0 


5 X 10-3 


0.120 


0.804 


0.994 


1.128 


0.45 


25.0 


5 X 10-3 


0.120 


0.858 


1.047 


1.182 


0.50 


16.0 5 X 10-^ 


0.115 


0.912 


1.098 


1.230 



experiment [191 || 1.230 | 1.270 | 1.170? 



TABLE 5. Results for the axial-vector (J^ ~ 1^^) mesons. 
All units are GeV. 



CO 


D 


mu 


ms II M„i7 


M„5 




0.40 


45.0 


5 X 10-3 


0.120 


0.917 


1.117 


1.253 


0.45 


25.0 


5 X 10-3 


0.120 


0.918 


1.124 


1.270 


0.50 


16.0 5 X 10-3 


0.115 


0.927 


1.140 


1.292 


ext^eriment [191 || 1.230 


1.270 


1.282 



ties are not too sensitive to the conjectural parameter dependence of the timelike quark 
propagator functions. Presumably meson properties whose computation involves larger 
timelike momenta will exhibit a stronger sensitivity. 

Already from table H] we observe that by increasing 0) the predicted mass of the 
J^'-' = meson with pion flavor quantum numbers approaches the empirical mass. We 
therefore further increased CO according to the rules discussed above. For O) ~ O.SGeV 
we reproduced the empirical value for the mass in that channel. However, this happened 
at the expense of significantly lowering fx and loosing the proper description of the 
vector mesons. We recall that the parameter o) has a physical interpretation as the 
location of the maximum of the interaction. Thus O) = O.SGeV seems intuitively too 
large for low-energy hadron physics and an unsatisfactorily description of the and 
1- mesons comes without surprise. 

For non-diagonal flavor structures such as us, charge conjugation actually is not a 
sensible quantum number and the corresponding axial vector mesons 1++ and 1+- may 
mix. In table 1^ we present the results obtained from the full calculation that combines 
the Dirac decompositions Q and dlOb . Since our Bethe-Salpeter formalism only yields 

TABLE 6. Results for the axial-vector mesons allowing for 
mixing of the Dirac structures in eqs. (|9j and (I10> . All units 
are GeV. 



CO 


D 


m„ 


ms II M„u 


M„5 




0.40 


45.0 


5 X 10-3 


0.120 


0.807 


0.990 


1.131 


0.45 


25.0 


5 X 10-3 


0.120 


0.861 


1.040 


1.185 


0.50 


16.0 


5 X 10-3 


0.115 II 0.915 


1.085 


1.233 



the lowest mass eigenstate within a given channel, the results presented in table|6l should 
be compared to those in table IH The tiny changes for the flavor diagonal mesons are 



— pseudoscalar 

— scalar 

— vector 

— axialvector (1 +-) 
axialvector (1 ++) 



Qt ^ \ , \ ^ 1 

0.4 0.8 

FIGURE 2. Meson masses as a function of the (equal) quark mass, mo. CO = 0.5GeV,D = 16GeV 
All units are GeV. 



TABLE 7. Predicted masses of cc-meson states. © = 0.4, D — 45.0, — 
1.125. nic is fitted approximately from the Tj^ mass. 



II 0- 


1- 


0+ 


1++ 


1+- 


Mcc II 2.97 


3.13 


3.32 


3.38 


3.31 


experiment || rjc: 2.98 


J/xj/: 3.10 


Xco: 3.42 


Xci-- 3.51 


? 



numerical artifacts. Surprisingly the changes for the non-diagonal flavor structure are 
also only of the order 1%. This suggests only a small mixing between the 1++ and 
1+^ states with the flavor structure us and the 1++ and 1^ channels represent good 
approximations to the actual eigenstates. 

We can extend our model beyond the light flavors up, down and strange. The only 
modification is the increase of the current quark mass, mq. In Fig.|2lwe show how the 
meson masses increase as the (equal) quark masses are increased into the charm sector 
Me — 1.125GeV. Exactly the same numerical code is used to construct these solutions 
to the Schwinger-Dyson and Bethe-Salpeter equations as for the light flavors. Clearly 
seen is the smooth way the masses increase from the chiral limit (mg = 0) into the 
heavy quark sector (mo = rric). This represents a convincing indicator for the stability 
of our technique. The cc-meson masses can be loosely extracted (table |7l) and the data 
are surprisingly well reproduced. The lack of the correct UV behavior for the gluon 
is seemingly at odds with the scales present. However, the present results suggest that 
the Bethe-Salpeter equation is capable of describing all the angular momentum states 
equally well in the charm quark sector. 

SUMMARY AND OUTLOOK 

In this talk we have presented a study of the low-lying mesons as quark-antiquark 
bound states in a covariant approach using an effective interaction. This interaction is 




characterized by gluon exchange with the gluon propagator being dressed by a GauBian 
shape function. The interaction is completed by the quark-gluon vertex that we take to 
be the tree-level perturbative one. In this manner the rainbow-ladder approximation 
to the system of Schwinger-Dyson and Bethe-Salpeter equation accounts for chiral 
symmetry. With this effective interaction, we have then consistently treated this system 
of integral equations by precisely implementing the quark propagator functions that 
solve the Schwinger-Dyson equations into the Bethe-Salpeter equations. Once the 
effective interaction exceeds a certain strength, the Schwinger-Dyson equations exhibit 
dynamical chiral symmetry breaking and the pseudoscalar mesons emerge as would-be 
Goldstone bosons. We have then used observed properties of the pseudoscalar mesons to 
determine the model parameters. The kaon decay constant represents a model prediction. 
It turned out to be in good agreement with the empirical data. Furthermore our results for 
the vector meson masses match the experimental data. The situation in the scalar channel 
is less satisfying. As we solely consider the mesons as bound states of quark-antiquark 
pairs, it is not surprising that the mass eigenvalues increase with the strangeness content. 
On the other hand it is astonishing that for current quark masses, mq > 0.2GeV, the 
lightest scalar mesons turn out to be heavier than the lightest vector mesons. When 
discussing these results it must be noted that the role and structure of the scalar mesons 
is still under intense debate. In particular, the question whether they should indeed be 
considered as quark-antiquark bound states is not yet completely resolved. There are 
indications, see e.g. ref. [1] and references therein, that the scalar meson masses should 
actually decrease with the strangeness content of these mesons. This can be understood if 
these mesons are considered as 2quark-2antiquark bound states in the sense of diquark- 
antidiquark systems 

As an outlook we mention that there is an elegant way to extend the present model 
to incorporate such degrees of freedom. The Bethe-Salpeter treatment can be straight- 
forwardly extended to study bound states of diquark-antidiquark pairs, once a binding 
mechanism is established. This could either be achieved by a gluon exchange similar to 
eq. dl]) or by quark exchange between a quark and a diquark. The latter approach has 
been intensively studied and the corresponding vertex is known from modeling baryon 
properties [20]. It will also be interesting to see whether these additional degrees of free- 
dom will also affect the mass predictions for the axialvector mesons that currently tend 
to be on the low side. Investigations in this direction are in progress. 
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